Role of finite layer thickness in spin-polarization of GaAs 2D electrons in strong 

parallel magnetic fields 
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We report measurements and calculations of the spin-polarization, induced by a parallel magnetic 
field, of interacting, dilute, two-dimensional electron systems confined to GaAs/AlGaAs heterostruc- 
tures. The results reveal the crucial role the non-zero electron layer thickness plays: it causes a 
deformation of the energy surface in the presence of a parallel field, leading to enhanced values for 
the effective mass and g-factor and a non-linear spin-polarization with field. 
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PACS numbers: 73.50.-h, 71.70.Ej, 73.43.Qt 



The spin-polarization of an interacting, dilute two- 
dimensional (2D) carrier system has been of interest 
for decades. It has long been expected that because 
of Coulomb interaction the product g*m*, which de- 
termines the spin susceptibility of the 2D system, in- 
creases as the 2D density (n) is lowered and eventually 
diverges as the system makes a transition to a ferromag- 
netic state at sufficiently low n Q, (g* and m* are 
the carrier Lande g-factor and effective mass, respec- 
tively). Recently, there has been much renewed inter- 
est in this problem, thanks to the availability of high- 
quality dilute 2D systems, and the belief that it may 
shed light on the controversial issue of a metal-insulator 
transition in 2D |3|. A technique commonly used to 
study the spin-polarization is to measure the response 
of the 2D system to a tilted or parallel magnetic field 
H H H 11 11 EjjlJEI The results of some of these 
measurements @,Sll3i however, appear to be at odds 
with what is theoretically expected 0, f° r a dilute, in- 
teracting 2D system that is otherwise ideal, i.e., has zero 
layer thickness and is disorder-free. In particular, when 
g*m* is deduced from parallel magnetic field at which the 
2D system becomes fully spin-polarized, then the experi- 
mental results for GaAs 2D electrons Q and holes 0, ITcj 
suggest a decreasing value of g*m* with decreasing n, 
opposite to the theoretical predictions. 



Here we report a combination of measurements and 
calculations for the parallel magnetic field-induced spin- 
polarization of 2D electrons at the GaAs/AlGaAs het- 
crojunction. The results highlight the importance of the 
finite thickness of the electron layer and the resulting de- 
formation of the energy surface .E(k||), where kp is the 
in-plane wave vector, that occurs in the presence of a 
strong parallel field. This deformation induces an en- 
hancement of both m* and g* and leads to a non-linear 
spin-polarization in a parallel field. We find that, once 
the effect of the finite layer thickness and interaction is 
taken into account, there is reasonable agreement be- 
tween the experimental data and calculations. 

We used five samples from three different wafers (A, 



B, and C). The samples were all modulation-doped 
GaAs/AlGaAs hetero junctions with n in the range 0.8 
to 6.5 x 10 10 cm -2 . Their low-temperature mobility var- 
ied depending on the sample and n; at n = 2 x 10 cm" 2 , 
it ranged from about 2 x 10 5 to 2 x 10 6 cm 2 /Vs. Sam- 
ples were patterned in either van der Pauw or Hall bar 
shapes, and were fitted with back- or front-gates. To tune 
n in samples A, Bl, and B2, following illumination with 
a red LED, we used front-gate bias; for B3 (in the range 
n < 4 x 10 10 cm -2 ) and C we used back-gate bias and no 
illumination. For B3, the highest density (n = 6.5 x 10 10 
cm~ 2 ) was obtained after illumination, followed by back- 
gating to reduce n to 4.5 x 10 10 cm -2 . Measurements 
were done down to a temperature of 30 mK, and a ro- 
tating platform was used to tune the angle between the 
applied magnetic field and the sample plane. 

Figure 1 summarizes our data taken on the different 
samples. Plotted are the values of g*m* j g&mf,, deter- 
mined from the parallel magnetic field, Bp, at which the 
2D system becomes fully spin polarized (m& = 0.067mo 
and <?b = —0.44 are the band effective mass and Lande 
g-factor for GaAs electrons; mo is the free electron mass). 
The parallel magnetic field (Bn) leads to the formation 
of two energy subbands, one for each spin, and sepa- 
rated by the Zeeman energy, Ez — |<7*|//b-Bm, where \xb 
is the Bohr magneton. The 2D system becomes fully 
spin-polarized above a field Bp at which Ez equals the 
Fermi energy. The equality leads to an expression for Bp: 
Bp = (h 2 /2tthb) • ( n /\g*\ m *), from which we determine 
g*m* that are plotted in Fig. 1. 

The procedures we have used to experimentally de- 
termine Bp have been described elsewhere 0, here 
we give a brief summary. We determine Bp from two 
independent sets of experiments: Shubnikov - de Haas 
(SdH) measurements in a nearly parallel magnetic field, 
and magnetoresistance measurements in a strictly par- 
allel field. In the first type of experiment, we apply a 
constant magnetic field (B to t) whose initial direction is 
parallel to the 2D electron plane, and then slowly rotate 
the sample while recording the sample resistance as a 
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FIG. 1: Values of g*m* / gtrrib, determined from the paral- 
lel magnetic field Bp at which the 2D electrons in GaAs 
become fully spin-polarized, are shown as a function of 2D 
electron density. Results are shown for five samples with low- 
temperature mobilities (at n — 2 x 10 10 cm -2 ) of A: 4, Bl: 
16, B2: 16, B3: 7, C: 2 x 10 5 cm 2 /Vs. The closed symbols 
represent Bp determined from SdH measurements in a nearly 
parallel magnetic field while open symbols are for Bp from 
magnetoresistance measurements in a parallel field. 



function of the angle between the 2D plane and the field 
direction. If we limit ourselves to small angles, the field's 
parallel component (Bu ) remains almost constant (equal 
to B to t typically to better than 1%) during the rotation, 
while its perpendicular component (B±) changes suffi- 
ciently to probe the SdH oscillations. We then Fourier 
analyze the SdH oscillations to obtain the populations of 
the two spin subbands. These densities provide a direct 
measure of the spin-polarization of 2D electron system 
and allow us to determine the field Bp above which the 
system becomes fully spin-polarized. In the second type 
of experiment we measure the sample resistance as a func- 
tion of a magnetic field applied strictly in the 2D plane. 
As shown elsewhere 0,0, the in-plane magnetoresistance 
shows a marked change in its functional form at the field 
Bp. For a given sample, the value of Bp obtained from 
the two types of experiments (open and closed symbols 
in Fig. 1) are in agreement. 

Data of Fig. 1 illustrate that the product g*m*, de- 
duced from Bp as described in the last two paragraphs, 
deviates substantially from what is expected for an ideal, 
interacting 2D electron system, namely a monotonically 
increasing g*m* as n is lowered The data also re- 

veal that the measured g*m* is sample dependent and 
not a unique function of n. A possible reason for this 
non-uniqueness may be the sample disorder that indeed 
varies between different samples. An examination of the 
data, however, argues against this hypothesis: consider- 
ing the data at a given density, it is clear that there is no 
simple trend linking the sample disorder, as deduced from 
the low-temperature mobility, to the measured g*m* . As 
we demonstrate below, another factor that renders the 



experimental 2D electrons non-ideal, namely their finite 
layer thickness, appears to be responsible for the sample- 
dependent g*m* and the difference between the observed 
and expected density dependence of g*m*. 

In the presence of a large Bu , when the magnetic length 
(= y/h/eB\\) becomes comparable to or smaller than the 
thickness of the electron layer, the energy surface B(kii) 
of the electrons gets deformed in the in-plane direction 
perpendicular to B\\ . The deformation leads to an in- 
crease of the in-plane effective mass, to*, which, in second 
order perturbation theory, is given by |12( 
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where the sum runs over all excited subbands j and z 
is the quantization axis jl_3| . Data of Fig. 2 provide an 
experimental demonstration of this effect in our samples. 
In Fig. 2(a) we show, as a function of Bn, the measured 
magnetoresistance of samples A and Bl, both at a den- 
sity of 2.7 x 10 10 cm' 2 . The magnetoresistance for each 
sample shows a clear change in its dependence on Bn at 
a field marked by a vertical arrow as Bp. As demon- 
strated previously 0, Q , the field Bp marks the onset of 
full spin-polarization. Note in Fig. 2(a) that Bp is larger 
for sample Bl than A even though they have the same 
density. We also measured to* in the two samples as a 
function of Bn , as shown in Fig. 2(b) . We determined 
to* from the temperature dependence of the amplitude of 
the SdH oscillations, measured as the sample was slowly 
rotated in an almost parallel field. We performed a stan- 
dard analysis, fitting the amplitude of the SdH oscilla- 
tions (Ai?) to the Dingle formula, AR ~ £/ sinh£, where 
£ = 2ir 2 kpT /Tllli c and u> c = eB/m* . It is clear in Fig. 2(b) 
that to* for both samples exhibits a strong enhancement 
with increasing Bit, consistent with Eq. JTJ. Moreover, 
to* has a larger enhancement for sample A than for Bl. 
This stronger enhancement correlates with the smaller 
Bp measured for sample A. We believe that the main 
difference between the two samples in Figs. 2(a) and (b) 
is that sample A has a larger layer thickness: because 
of a decrease in subband separation, the larger thickness 
leads to the larger enhancement of to*, consistent with 
the smaller Bp for sample A compared to Bl. 

To further substantiate the connection between layer 
thickness and to* enhancement, in Fig. 2(c) we show data 
for sample A at a higher carrier density of 4.2 x 10 10 
cm~ 2 . The measured to* enhancement is smaller for the 
higher n state. This is consistent with layer thickness 
being responsible for the to* enhancement: as we use a 
more positive front-gate bias to increase n in this sample, 
the electron wavefunction is squeezed more towards the 
interface so that the layer thickness is reduced [see the 
insets to Fig. 2(c)], consistent with the smaller measured 
m* enhancement. 

In order to quantitatively understand the experimental 
data, we have done self-consistent density-functional cal- 
culations of the subband structure in the presence of Bu . 
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FIG. 2: (a) Parallel-field magnetoresistance of samples A and 
Bl at 2.7 x 10 10 cm" 2 , (b) Effective masses measured for the 
two samples of panel (a), (c) Effective masses measured for 
sample A at two different densities. The density is tuned via a 
front-gate bias which leads to a narrowing of the wavefunction 
at higher density, as shown schematically in the insets. In 
all three panels, the fields Bp above which the 2D electrons 
become fully spin-polarized are marked by vertical arrows. 

We used the recent parameterization of the exchange- 
correlation energy by Attaccalite et al. Q- The en- 
ergy and length scales for electrons in a semiconduc- 
tor are characterized by the effective Rydberg and the 
Bohr radius according to the effective mass m* and the 
dielectric constant of the material. In our calculations 
it was crucial that m* was determined as a function of 
B\\ from the self-consistently calculated subband disper- 
sion E(k.u, Bu). These calculations confirm the qualita- 
tive trends expected from Eq. (JIJ . Band structure effects 
beyond the effective-mass approximation are not consis- 
dered here. We have checked that these effects are of 
minor importance In the calculations, the field Bp 
is defined as the smallest value of B» for which the fully 
spin-polarized configuration has the lowest total energy. 




FIG. 3: Results of calculations are shown for (a) spin subband 
densities n-f and nj,, (b) effective mass, and the g-factor for 
a GaAs 2D electron system (density 4.2 x 10 10 cm -2 ) with 
realistic finite layer thickness. Also plotted (circles) are the 
experimentally measured m* , n\ and ny at this density. 



Figure 3 provides an example of the results of the cal- 
culations for n = 4.2 x 10 10 cm~ 2 ; shown are (a) the 
spin-subband densities n-\ and rij, (b) m*, and (c) g*, 
as a function of B\i . The calculations were done using 
the parameters of sample A (spacer thickness and bar- 
rier height), assuming a p-type background doping of 
2.7 x 10 13 cm -3 , and a binding energy of 90 meV for the 
dopants (Si) in the barrier. These are reasonable values, 
consistent with our estimate of the unintentional (resid- 
ual) doping in our molecular beam epitaxy system and 
the binding energies quoted in the literature |l6l |. The 
calculations predict a non-linear but smooth increase of 
the spin-polarization as a function of Sii . 

In Figs. 3 (a) and (b) we have also included our mea- 
sured spin-subband densities and to*, determined from 
SdH oscillations in a nearly parallel magnetic field. Over- 
all, there is good agreement between the experimental 
data and calculations for both spin-polarization and m* 
as a function of Bu. The calculated Bp = 14.5 T agrees 
well with the measured Bp = 16 T, and is much smaller 
than Bp = 48 T expected for an interacting GaAs 2D 
electron system with zero layer thickness at n = 4.2 x 10 10 
cm~ 2 (see dashed curve in Fig. 4). 

To understand the density dependence of Bp for a 
given sample, we also calculated Bp as a function of n, 
which is tuned either with a front- or back-gate bias. 
For these calculations, we kept the sample parameters 
fixed, and only changed the boundary conditions for the 
Hartree potential thus simulating the effect of the gate 
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FIG. 4: Comparison between the measured Bp for samples 
A and Bl (same symbols as in Fig. 1) and calculations that 
take the finite layer thickness of the GaAs 2D electrons into 
account (solid curve) . The density was varied via a front-gate 
bias in both the experiments and calculations. For contrast, 
we also show Bp from two sets of calculations that are based 
on Ref. and are done for interacting, zero layer thickness 
2D electrons in GaAs (see text). The lower inset shows the 
calculated curves of the main panel at lower densities. The 
upper inset shows g*m* based on different calculations. 

bias. The results for the case where the front-gate is used 
are shown in Fig. 4 (solid curve), along with the experi- 
mental data for samples A and Bl. There is reasonable 
agreement between the calculations and the experimental 
data 

In Fig. 4 we also plot Bp vs. n (dashed curve) calcu- 
lated based on Ref. 0] where an interacting 2D electron 
system with zero layer thickness is assumed, as well as 
" Bp" (dotted curve) determined from the calculated spin 
susceptibility in the limit B\\ = and assuming a lin- 
ear spin-polarization as a function B\\ . Three noteworthy 
trends are observed in Fig. 4. First, for n > 0.5 x 10 10 
cm~ 2 the dashed and dotted curves are very close to each 
other, meaning that the spin polarization of an ideal, 



zero-thickness 2D system is approximately (within 5%) 
linear with B^ in this density range or, equivalently, the 
Zeeman splitting can be expressed well in terms of an ef- 
fective g-factor independent of Bu . This implies that the 
finite layer thickness is the key factor that leads to the 
observed non-linearity in spin-polarization with Bu and 
the resulting reduction of Bp . The mechanism responsi- 
ble for this non-linearity can be summarized as follows: 
induces an increase of m* due to the finite layer thick- 
ness. This has a twofold effect. It directly reduces Bp 
because Bp oc 1/m*. Furthermore, the increase of m* 
reduces the effective Bohr radius and thus increases r s , 
the average electron spacing measured in units of effec- 
tive Bohr radius. The increase in r s in turn yields an 
increase of g* due to the Coulomb interaction. Second, 
the solid and dashed curves merge as the density is low- 
ered, consistent with the expectation that, because of the 
smaller Bp, the finite layer thickness induced g*m* en- 
hancement becomes less important. Third, at ultra- low 
densities the dashed and dotted curves start to diverge 
(see the lower inset), implying that the interaction alone 
can induce a significant non-linearity of spin-polarization 
with B\\ in a very dilute 2D system 0. 

Finally, if we use the expression Bp = (h 2 /2tt^lb) • 
(n/\g*\m*) to determine g*m* from Bp as frequently 
done in the literature, we obtain the curves shown in 
Fig. 4 upper inset. These plots emphasize that "g*m*" 
deduced from Bp for a 2D system with finite layer thick- 
ness (solid curve) is significantly enhanced with respect 
to the ideal 2D system and can show a non-monotonic 
dependence on n. The results therefore caution against 
extracting values for g*m* in the limit of zero magnetic 
fiel d from measurements of Bp at large parallel fields 
[18j. Moreover, Bp and g*m* are not unique functions 
of n; they depend on the electron layer thickness which 
in turn depends on sample parameters and experimental 
conditions. 
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